Thermal transport in the Fahcov-Kimball model 
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We prove the Jonson-Mahan theorem for the thermopower of the Fahcov-Kimball model by solving 
explicitly for the correlation functions in the large dimensional limit. We prove a similar result for 
the thermal conductivity. We separate the results for thermal transport into the pieces of the heat 
current that arise from the kinetic energy and those that arise from the potential energy. Our method 
of proof is specific to the Falicov-Kimball model, but illustrates the near cancellations between the 
kinetic-energy and potential-energy pieces of the heat current implied by the Jonson-Mahan theorem. 
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I. INTRODUCTION 



The Jonson-Mahama theorem shows that there is a 
simple relation between the transport coefRcient for the 
electrical conductivity and that needed for the ther- 
mopower. The relation is that the integral for the L12 
coefficient has an extra power of frequency in the in- 
tegrand than the Ln coefficient. This result has been 
known for many yeara3 for a noninteracting system — 
the Jonson-Mahan theorem generalizes this result for a 
wide class of many-body systems (including the Falicov- 
Kimball model, the static Holstein model, and the pe- 
riodic Anderson model, but not including the Hubbard 
model or quantum Holstein model). 

We use the exact solution of the Falicov-Kimball model 
in the large-dimensional limit to provide an alternate 
derivation of the Jonson-Mahan theorem by explicitly 
evaluating all relevant correlation functions needed for 
the thermal transport. Our exact analysis also allows us 
to separate the contributions to thermal transport that 
arise from the kinetic energy and the potential energy 
pieces of the heat current. These results provide an inter- 
esting interpretation of thermal transport in correlated 
systems. 

In Section H we develop the formalism for deriving the 
dc conductivity, the thermopower and the thermal con- 
ductivity. We derive exact results for the relevant corre- 
lation functions and use them to prove the Jonson-Mahan 
theorem and its generalization to the thermal conductiv- 
ity. In Section HI we provide numerical results for the 
thermal transport illustrating the different contributions 
to the thermal coefficients for a number of illustrative 
cases. Conclusions are presented in Section IV. 



II. FORMALISM FOR THE THERMAL 
TRANSPORT 



The Hamiltonian foj: our system is the spin-onc-half 
Falicov-Kimball modelQ 



H = 



2Vd 



E 



(1) 



where (cia) is the electron creation (annihilation) op- 
erator for an electron at site i with spin a, Ef is the 
energy level of the localized electrons, Wi is a variable 
that equals zero or one and corresponds to the localized 
electron number, and U is the interaction strength. The 
hopping integral is scaled with the spatial dimension d so 
as to have a finite result in the limitu d —> 00; we measure 
all energies in units of i* = 1. We work on a hypercu- 
bic lattice where the noninteracting density of states is a 
Gaussian p(e) = exp(— e~)/Y^. 

The Falicov-Kimball model can be saUfed exactly by 
employing dynamical mean field theoryO'El. Because the 
self energy T,{z) is local, the local Green's function satis- 
fies 



G(z) = / dep{e) 



1 



z + ^ — 'E{z) 



(2) 



with z anywhere in the complex plane (we suppress the 
spin index here). The self energy, local Green's function, 
and effective medium Go are related by 



Go\z)~G-\z) = ^iz), 



and the Green's function also satisfies 



G{z) = (1 - wi)Gq{z) + wi 



1 



G^'{z)-U 



(3) 



(4) 



Here wi is the average concentration of localized elec- 
trons, 

wl = 2eM-f3{Ef-^i)]z^{^l-U)Zl{^i-U)/z, (5) 

with Z = Z^{|l)Z^{^i) + 2exp[-/3(£;/ - m)]^t(m - 
U)Zi{^i-U) and 



ZA^^)^2ef'^/'Y[ 



iUJn + /-i - Aa(«W„) 



(6) 



1 



The factor of two arises from the spin degeneracy of the 
f-electrons and the constraint that no more than one f- 
electron is allowed on any site. The symbol Xa{z) is de- 
fined from the effective medium via X^{iu!n) — i^On + 
GQ]:{iujn), ujn — TrT{2n + 1) is the Fermionic Matsubara 
frequency, and (3 = l/T. The algorithm for determining 
the Green's function begins with the self energy set equal 
to zero. Then Eq. (||) is used to find the local Green's 
function. The effective medium is found from Eq. and 
the localized electron filling from Eq. (H). The new local 
Green's function is then found from Eq. (|) and the new 
self energy from Eq. (||). This algorithm is repeated until 
it converges. 

Transport properties are calculated within a Kubo- 
Greenwood formalism^. This relates the transport co- 
efficients to correlation functions of the corresponding 
transport current operators. We will deal with two cur- 
rent operators here — the particle currentEl 



J = 



(7) 



(where the velocity operator is Vg = Vqe{q) and the 
Fourier transform of the creation operator is = 

J2j exp[i(7 • Rj]c|/N) and the heat current!'! 



qa 
U 



W{q - q')[Vq + Vq,]cl^Cq'a, 



(8) 



[where W{q) = J2j cxp(— i(7-Rj)wj/A^]. The heat current 
can be broken into two pieces: (i) a kinetic-energy piece 
jg which is the first term in Eq. (H) and (ii) a potential 
energy piece jg which is the second term in Eq. (^. 

The particle current is defined by the commutator of 
the Hamiltonian with the polarization operatorJj When 
evaluated on a lattice with nearest-neighbor hopping, 
one finds factors that involve the weighted summation 
of the nearest-neighbor translation vectors S weighted by 
phase factors exp{iq ■ S), which yield the velocity opera- 
tor terms above. The definition of the heat-current op- 
erator is more involved, and requires the Hamiltonian 
to be separated into operators hi that involve the site 
i (in decomposing the kinetic- energy operator into "lo- 
calized" pieces, one symmetrically assigns half of the 
cJ^Cjo- -|- Cj^Ci„ term to site i and half to site j). These 
operators can be combined with the position operator 
to construct an "energy" polarization operator ^^R^/ii, 
which is commuted with the Hamiltonian to determine 
the energy current operator; the heat current operator is 
just this energy current operator shifted by the chemi- 
cal potential multiplied by the number current operator. 
Important operator relations between the heat-current 
operator and the particle-current operator are described 
fully below [see the discussion around Eq. (IgT 



The dc conductivity cr, thermopower S and electronic 
thermal conductivity k can all be determined from rele- 
vant correlation functions of the current operators. We 
define three transport coefficients Ln, L12 = L21, and 
L22. Then 



T 



S = 



jelTLn' 



and 



L 



L12L 



21 



22 



Lu 



(9) 



(10) 



(11) 



The transport coefficients are found from the analytic 
continuation of the relevant "polarization operators" at 
zero frequency 

2 — 

ill = lim Re— Lii(i^), 

L,^{iui)^t:T dTe'-^^Trjl{T)jp{Q)), (12) 
Jo 

where vi — 2ttTI is the Bosonic Matsubara frequency, 
the r-dependence of the operator is with respect to the 
full Hamiltonian in Eq. (p, and we must analytically 
continue Ln^ivi) to the real axis Lii(^) before taking 
the limit v —^ 0. Similar definitions hold for the other 
transport coefficients: 



and 



L12 

Li2{wi) 

L22 - 
£22 (if/) = 



lim Re— Li2(j^), 
Jo 



(13) 



lim Re-L22(i^), 

nT I dTe-'-(T.4jr)jQ^(0)). (14) 



In all of these equations, the subscripts a and /? denote 
the respective spatial index of the current vectors. 

We begin with a derivation that shows the analytic 
continuation for the conductivity. Substituting the def- 
inition of the particle current operator of Eq. (^ into 
Eq. (0) for Lii yields 



qq'aa' 

X {TrcliT)Cq,iT)cl,,,iO)Cq^,-m. 



(15) 



The correlation function can be determined from Dyson's 
equation which relates the dressed correlation function to 
the bare correlation function via the irreducible charge 



2 



vertex. Since the charge vertex is local in the infinite- 
dimensional limit, it is an even function of momentum, 
and any sum over momentum that is weighted by just 
one factor of will vanish. Hence, the dressed cor- 
relation piunction is equal to just the bare correlation 
fmictionEEl (note that the contractions of the operators 
at equal times also vanish when summed over momen- 
tum). This produces 

Lii{ivi) = -vrT / dre^'^'" Vv,aV,/5G,.(T)G,.(-r). 
Jo 



qa 



(16) 



Now, we introduce the Fourier transform of the Green's 
function G{t) = T exp(— iCi;„r)G„ with G„ = 



dr exp(ia;„r)G(T). Substituting into Eq. (16), allows 



us to perform the integral over r. This finally produces 

n qa 

(17) 



(b) 
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Freericks and Zlatic, Phys. Rev. B, Figure 1 (b) 

FIG. 1. Contours used in the analytic continuation: (a) 
contour needed for the Matsubara frequency summation and 
(b) deformed contour to lines parallel to the real axis. 
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Freericl^s and Zlatic, Phys Rev B, Figure 1 (a) 



The next step is to perform the analytic continua- 
tion from the imaginary to real axis. The procedure 
is standards. We first write the summation over Mat- 
subara frequencies as an integral over the contour G 
shown in Fig. 1(a) which has contributions at the poles 
of the Fermi function /(w) = 1/[1 + exp(/3a;)] which lie 
at the Fermionic Matsubara frequencies. The contours 
are then deformed to lines parallel to the real axis, with 
the Green's functions evaluated with either retarded (R) 
or advanced (A) functions. The result is 



Lllim) 



T 

2ijc 



T 



2i 



/OC 
-oo 



AG'^,(u:)-G%C^)\G%{u: + ^v^) 



q0 



rp PCX) 

-^-^ qa 

xGfA^~^,,l)[GfAoJ)~GfAuJ)]. (18) 

The analytic continuation is performed by first rewriting 
/(w — ii/i) = f{uj), then taking ivi ^ v + iS and shifting 
the integration variable w — > w + i' in the second integral. 
Then, using the definition for in, we finally arrive at 



Lll = lim -— / dujy^ VqaVq/S 



X Rc <^f{uj)Gq„{uj)Gqaiuj + 1^) 
^ f{u; + iy)G;A^)G*qA^ + ,.) 



3 



(19) 



Since Gqa-i^) = + — — e(g)], we can per- 

form the summation over q directly. Because e(g) is an 
even function of q and is odd, we must have a = f3. 
Converting the fraction into the integral of an exponen- 
tial, thes, allows the summation over q to be performed 
directlyllil. The summation over q can be written as an 
integral over energy with a weighting factor of p{e)t*'^ /d. 
This yields 



Lii = lim -1^— 7<5q/3 



dco 



dtp{e) 



X Re<yf {uj)Gqc{uj)Gqa{LO + v) 

- f{LO + u)G*q,{uj)G*qAuJ + u) 

- [J{uo) - fiio + iy)]Gl,iu;)Gq.iu; + (20) 

If we define (Tq = e^i*^/ {2d), and we perform the integral 
over e, we arrive at 



Lii = lim '^-^Sap [ duj 

V 

^ G*(i^)~G{uj + v) 
V + — S(a; + v) 



G'(^) - G{uj + v) 
V + S(a;) — Y^{u} + v) 



(21) 



The final step is to take the limit of ^ 0. Using the 
facts that 



lim 



V div ' 



(22) 



and 



G{u) - G{uj + v) 
i'^o V -\- S(a') — E(a; -|- v] 



lim 



= -2 + 2[tj + /i-I](tj)]G(t^), 



produces our final result for Ln 



ill 



Tctq 



do; 



dm 

doj 



(23) 



(24) 



with the relaxation time t{lu) defined by 

r{cj) = ^^^^ + 2 - 2Re{[c. + /i - n^)]G{uo)}. (25) 
lmL(a;) 

This result appears differept from that originally derived 
for the hypercubic latticeEJ, but Eqs. (|2J) and (25) do 
yield the same result; the form presented here has the in- 
tegral over the noninteracting density of states performed 
exactly. 

The derivation for the transport coefficient L12, needed 
for the thermopower, proceeds in a similar fashion. We 



can divide the heat current into two pieces, one cor- 
responding to the kinetic energy and one correspond- 
ing to the potential energy. This allows us to write 
L12 — L^2 + ^12- The derivation for the kinetic energy 
piece follows exactly like the derivation for Ln except 
there is an extra factor of e — /i that appears in Eq. (|2(j). 
The integral over e can then be performed straightfor- 
wardly producing 



= hm 



Tar 



Sal: 



duj 



X Re<^ - 



[lo - T.{uo)]G[uj) - [uj + V - Y.{uj + v)]G{lo + v) 
V -t- Yj(lo) — S(a; -f v) 
\uj - E*(w)]G*(u;) - [cj + J/ - E(u; v)\G{uj + v) 
i^ + E*(w)-i;(tj + j/) 

Evaluating the limit 1/ — > is simple. The final result is 

Tctq 



(26) 



^12 — 2 



'^"^^ Ht"^ " ReI](w)]T(w) 



2ImE(w)Im[(w + ^ - S](w))G(w)]}, 



(27) 



with t(lo) defined in Eq. (psj). 

The derivation of the potential energy piece is much 
more involved. The first step is to replace the 
momentum-dependent operator VF(q — (?') by its Fourier 
transform. Simplifying the expression for yields 



ttTU 



/3 



dre' 



qq'aa' j 



2 JO 

e-^'^'-^' {Trw,cl{r)cq, {r)cl,^, (0)c,,, (0)) 
^«'-^^(T.z«,ct (r)c,.(r)ct (0)c,v'(0)) . 



(28) 



Noting that the wj operator commutes with the 
Fermionic operators, allows us to use Wick's theorem to 
rewrite the terms in the square bracket as 



N ^ 



{TrW,cl{r)c,40))Gq4T) 



{TrWjCqAT)clmGq,i-T) 



Sqq'dcrcr', (29) 



where we have taken the appropriate contractions (note 
the velocity operators guarantee that we need not worry 
about any vertex corrections) . The correlation functions 
in Eq. ( P9[) can be evaluated by taking the derivative 
with respect to the components of an infinitesimal field 
— ^jhjWj. These correlation functions have a factor 
of exp[— /3i/] in the numerator and a factor of Z in the 
denominator. In addition, the r-dependence of the oper- 
ators arises from factors of exp[±TiJ]. Since the operator 
Wj commutes with all Fermionic operators, it is easy to 
verify that the expression in Eq. (E^) becomes 



-E 



4 



Gn 



T- 



Sqq'Saa', (30) 



and 



which follows by first removing the Wj operator through 
the derivative, then expressing the Fermionic operator at 
site J through a Fourier transform, and finally evaluat- 
ing the Fermionic averages. Substituting this result into 
Eq. (p8|) then yields 



7tTHJ_ 
2N 



{[T 



qa j 

(Wj)] Gqa {iUJn )}Gqa- {iUJn+l ) 

d 



+ Gqa{iuJn)[T-Qj^ + {Wj)]Gqa{iuJn+l) 



(31) 



The derivatives need to be computed. Writing the 
momentum-dependent Green's function as a Fourier 
transform 



Gqaiiuin) 



1 

TV 



(32) 



and using the identity 



Gycr(zCJ„) = y2 Gika{ii^7i)G^i]^{iuJn)Gljcr{iuJn), (33) 



kl 



(with G ^ the matrix inverse of G) allows us to compute 
the derivative as 



_d_ 



Gqa{iLOn) 



1 



N 



Gjja (iljJn)- 



(34) 



But in a homogeneous phase, the derivative of the local 
self energy with respect to the local field, and the local 
Green's function are both independent of the site j, so 
we finally arrive at 



dGqgjiuJn) 

dh 



91] 

Gna—g^Gqaiicdn)- (35) 



Since the self energy depends only on G„ and wi, the 
derivative can be computed by taking partial derivatives 
and using the chain rule 



dT,na- dwi 

dioi dh 



dh 



1 — GJ 



(36) 



no- dG„ 



Each 
direct! 



the derivatives in Eq. ( pq ) can be found 



U 



dwi 1 + G„,(2E„, -[/)' 
dwi wi {1 — wi) 



dh 



T 



(37) 
(38) 



1 — G,, 



(1 + Gruj^na )(l + G„,[I]„,-?7]) 



'dG„ 



1 + G„,(2S] 

ri 



u) 



(39) 



Substituting these derivatives into Eq. ( p6| ) and perform- 
ing some straightforward simplifications that involve the 
quadratic equation that the self energy satisfiesQ finally 
yields 

[T^ + {w)]Gq„{iu^) = ^Gq„{lU^). (40) 

Now we are ready to perform the analytic continuation. 
First we substitute Eq. into Eq. (|l]) and we note 
that the sum over j cancels the factor of 1 /N 



n qcr 



X [T^aiiLUn) +^aiiLUn+l)]Gqa{it^n)Gqaiit^n+l)- (41) 

Next, we rewrite the sum over Matsubara frequencies as 
a contour integral and perform the analytic continuation 
in the exact same way as before. If we then evaluate 
we find 



L^2 ^ lim 



duj 



dep{t) 



X Re /(cj){[S(cj) + S(cj + u)]Gq{uj)Gq{uj H- u) 

- [E*(lj) + I](l^ + v)]Gl{uj)Gq{Lu + v)} 

+ f{u + v){[Y.*{uj) + E(w + v)]Gl{uj)Gq{Lo -f v) 

- + i:*{oj + v)]G:{u)G*^{u: + v)} 



(42) 



Now the integral over e can be performed and the limit 
V ^ can be taken. It becomes 



-'^12 - -jT 



duj — 



dm 

duj 



RcT.{uj)t{uj) 



(43) 



+ 2ImE(w)Im[{cj + ^ - S(tj)}G(tj)] 

and ( ^ ) yields the Jonson 
df {oj) 



Adding together Eqs. (| 
Mahan result of 



L12 — 



duj — 



duj 



t{uj)uj. 



(44) 



Our final derivation is for the thermal conductivity co- 
efficient L22. Like before, we separate this into pieces 
corresponding to the kinetic energy and the potential en- 
ergy: L22 = L^2^ + L^/' + L^2^ + L^^. Due to the sym- 
metry of the terms, we have = ^22^- The kinetic 
energy piece is simple to calculate. Like in our derivation 
for L{^, the steps are identical to the derivation for Ln 
except we have an extra factor of (e — /.t)^ in Eq. (pT 



5 



Performing the integration over e and collecting terms 
finally yields 



tKK _ Tcfo 
^22 - — 



duj — 



df{Lo) 



The derivation for = ^22^ is identical to that of 

except we have an extra factor of (e — /i) in Eq. (42). 
Performing the integration over e then produces 

X |[cj - Rci;(Lj)]Rci;(cj)r(tj) - IniG(cj)ImE(cj) 
+ 2[ImE(w)]2 - 2[ImE(w)]2Re{[w + ^ - I](w)]G(w)} 
+ 2[uj - 2ReS](w)]ImI](w)Im{[w + /i - S](w)]G(w)} 



The final term we must evaluate is i^2^. This is the 
most complicated term to evaluate and we are unable to 
do so following the same strategy as employed in the 
derivation. Instead, we proceed by an alternate method 
based on the equation of motion (EOM) technique. The 
EOM for the Fermionic creation and annihilation opera- 
tors (in the momentum basis) are 

d 



-4(r) = [e(g)-/i]ct,(r) + C/^T4^(fc)4^,,(r), 



(47) 



and 

d_ 
57 



Cga(-r) = -{t{q) - [l\Cqcj{T) - [7 ^ W(fc)Cq_fe^ (t). 



(48) 

These EOMs can be employed to express the correlation 
function of the heat-current operators in terms of deriva- 
tives with respect to imaginary time as shown below 



tPP 
^22 



{ivi) 



dre' 



E 



qq' q" q'" aa' 
i'Vqa + Vg///Q)(v,//3 + ^q"p) 

X {T,W{q - q"')W{q' - q")c\Ar)cq-'' ^r) 
4,,,(0)c,..,(0)) 



lim 



ttT / dre"-'^ V Vg,v,,^ lim li 



X [{1(5.. - - (6, - n)}cl,{r")c,,,iT'")]). (49) 



Now each of the operator averages can be expressed in 
terms of Green's functions, since the velocity factors 
guarantee there will be no vertex corrections. Noting 
further, that the integrals will only contribute if a = /3 
finally yields 



X |[w-ReS](w)]V(tj)+ImG(tj)I 

- 2[ImE(w)]2 + 2[ImI](w)]2Re{[tj + /i - S(w)]G(w)} 

- 4[uj - ReS(w)]ImE(w)Im{[w -f /i - E(u;)]G(w)} 



•'0 



qa 

-^G,Ar)d^rGqA~r) 

+ (e, - ^i){G,^{T)^rG,,i-T) - drG,,{T)G,„i-T)} 
-{eq - llfGqa{T)Gq„{~T) 



(50) 



We need to be able to produce expressions for the deriva- 
tives of the Green's functions. We do so by first writing 
the Green's function as a Fourier series over the Matsub- 
ara frequencies, and then taking the derivative into the 
^^^tsubara summation. This may appear to be math- 
ematically unsound, but we do so by adding and sub- 
tracting the known 5(t) behavior of the derivative, to 
regularize the summation. Since we are interested only 
in < T < /9, this procedure has no convergence issues. 
Likewise, one is also able to take the second derivative in 
this fashion. We find for < r < /3 

drGqair) = - {eq - fl)Gqcr{T) 



— lUJnrT 



dlGq„{T) ^ +{eq- ^ifGq,{T) 

+ (e,-M)r^e--' 



iujjn + M - 5]„ 



(51) 



with similar formulae for Gqa{—T). Substituting the 
derivatives from Eq. (51) into Eq. (|50|), and then sim- 
plifying the result finally produces 

^r2^(-o = "^EE^^. 

n qa 

{T.„{iuJn) + 'il„{iuj.n+l)YGq„{iuJn)Gqa{iuJn+l) 



+ Y.a{iuJn+l)Gqa{iuJn) + T.a{iuJn)G qa{iuJn+l) 



(52) 



It is easy to understand the first terms in this expres- 
sion, as they are what one would naively recover when 
following the same Wick analysis that was done previ- 
ously for L^2- We have not been able to discover a direct 
operator-based derivation of the second terms, but they 
are critical for providing the right answer for L^^. Per- 
forming the analytic continuation and simplifying yields 
our final result 



6 



rPP 
^22 



duj 



duj 



X |[ReI](tj)]V(w) - 2[ImE(tj)]2 + ImI](w)ImG(w) 
+ 2[ImS(w)]2Re{[w + /i - T,{uj)]G{uj)} 
+ 4ReI](w)ImI](cj)Im{[cj + ^ - S(tj)]G(tj)} 



Holstein model, and the periodic Anderson modeB, but 
it does not hold for either the Hubbard model or the 
quantum Holstein model, a fact which does not appear 
to be widely known. The "polarization operators" then 
become 



(53) 



Summing together Eq. (|4^), twice Eq. ([46|), and Eq. (|5. 
gives the Mott form 



Tao 



duj 



duj 



(54) 



We can also generalize the original Jonson-Mahan ar- 
gument to prove relations between L21 and in and be- 
tween L22 and Li2. Our method is different from their 
proof and relies on the infinite-dimensional limit, but one 
could proceed in their fashion if desired. We begin with 
the generalized two-particle correlation function 



E v,„v,,^(T.ct,(r)c,.(r')4,,,(r")c,v'(r"')). (55) 

qq' aa' 

In the infinite-dimensional limit, the two-particle correla- 
tion function is expressed by just its bare bubble because 
the irreducible charge vertex has a different symmetry 
than Vq . Hence, we immediately learn that 



= -Y.^lJo.pG,,{T"' - r)G,.(r' - r"). 

qa 

But 

Gq^T)^ J dcuA{k,u;)e-'^^[l-f{uj)], 
for r > and 



Gqair) = f dujA{k,uj)e 



(56) 



(57) 



(58) 



for r < 0. Substituting into Eq. (56), then yields 



771 / I n m\ 



'2d 



dep{e) dui duj' 



X A{e,Lu)A{e,Lu )e' 



f{Lo)[l~f{u')]. (59) 



Using this function we can construct the relevant "polar- 
ization operators" . Recalling the EOM in Eqs. ( p7[ ) and 
(Eq) shows that 



lim 



1 



a7 



i]Y.^A^^^M^')=h(j). (60) 



The Jonson-Mahan theorem will hold for any Hamilto- 
nian that satisfies Eq. ( |60| ) (which follows from the rele- 
vant commutators and equations of motion). This iden- 
tity holds true for the Falicov-Kimball model, the static 



Zii=^r/ e^'''^F(T,T-,0,0), 



(61) 



for the conductivity. 



L12 = ttT 



1 



d 



d 



2 b7-a7^)^(^'^'«'0)' 



(in the limit where t' ^ t ) for the thermopower, and 



£2 



d_ 



d 



d 



d 



(63) 



(in the limit where r' ^ t~, t'" t"^ , and r" ^ 0+) 
for the thermal conductivity. Because of Eq. (|59|), the 
analytic continuation is trivial (one first converts from 
imaginary time to Matsubara frequencies and then per- 
forms the Wick rotation to the real frequency axis), and 
if we note the identity 

fiio) - ficu + i.) = -fiLu)[l - f{LU + v)][e-^'' - 1], (64) 
then we can easily compute that 



Lii = / dep{€) I duj 



for the conductivity, 

Li2 = [ d€p{e) I duj [- 



dio 



dm 

dhJ 



for the thermopower, and 
L22 = / dep{e) [ duj 



dm 

duj 



A^ie,Lo), (65) 



A\e,Lu)Lo, (66) 



A\e,Lu)io^ (67) 



for the thermal conductivity. This proves Mott's form for 
the thermal transport [since / dep{e)A'^{e,uj) = t{^)]- 



III. NUMERICAL RESULTS 

We present results first for a case where the filling of 
the localized electrons is a constant. We choose the sym- 
metric case of (w) = 1/2 and half filling pe = 1 for the 
electrons. We perform calculations for two cases: (i) 
U = 1 which is a moderately correlated metal and (ii) 
U — 2 which is a strongly correlated insulator. 



7 




CD 



CO 
CD 

o 

CL 

o 
E 

CD 



0.5 1 1.5 2 

Temperature T [t*] 

Freericks, and Zlatic, Phys. Rev. B, Figure 2 

FIG. 2. Thermopower for the case U = 1, wi = 0.5, and 
Pd = 1- The main panel shows the two different contributions 

from the kinetic and potential energy pieces of the heat cur- 
rent and the inset shows the net thermopower. Note how the 
two pieces axe large and nearly cancel to produce S, and how 
there is a sign change near T « 1.4. 
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Freericks, and Zlatic, Phys. Rev. B, Figure 3 

FIG. 3. Thermal conductivity for the case U = 1, wi — 0.5, 
and Pd = 1- The plot shows the different contributions 
from the kinetic and potential energy pieces of the heat cur- 
rent. Note how the thermal conductivity is essentially de- 
scribed by the kinetic-energy-only piece for moderate corre- 
lation strength. 




to produce the net thermopower (note the three order 
of magnitude difference in the scales for the main figure 
and the inset). The thermal conductivity also appears as 
expected. We can see that while the contributions from 
the potential energy are critical in determining the right 
thermopower, they have a relatively mild effect in the 
thermal conductivity for a moderately correlated metal 
(note how close the total thermal conductivity is to the 
kinetic-encrgy-only contribution) . 

As we increase the correlation strength, so that the in- 
teracting density of states has a gap and the system is a 
correlated insulator, the behavior of the thermal trans- 
port changes. 




0.5 1 1.5 2 

Temperature T [t*] 

Freericks, and Zlatic, Phys. Rev. B, Figure 4 

FIG. 4. Thermopower for the case U = 2, wi = 0.5, and 
Pd = 1- The main panel shows the two different contribu- 
tions from the kinetic and potential energy pieces of the heat 
current and the inset shows the net thermopower. Note how 
the two pieces axe large and neaxly cancel to produce S and 
how the thermopower appears to diverge as T becomes small 
(calculations run into numerical problems as T — » 0). 



The thermopower behaves as expected — it vanishes for 
large and small temperature, it has an electron-like peak 
at lower temperatures and a sign change at a temperature 
on the order of half the bandwidth. What is surprising 
is that there is such a large compensation between the 
kinetic and potential energy pieces of the thermopower 
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Temperature T [t*] 

Freericks, and Zlotic, Phys. Rev. B, Figure 5 

FIG. 5. Thermal conductivity for the case U = 2, wi = 0.5, 
and /9d = 1. The plot shows the different contributions from 
the kinetic and potential energy pieces of the heat current. 
Note how the potential energy terms become increasingly im- 
portant. 

The thermopowcr has the characteristic insulating be- 
havior here, with what appears to be a divergence as 
T — > 0. The divergence arises from the presence of a 
gap in the single-particle spectrum — both Ln and L12 
approach zero exponentially in T (with the same expo- 
nent), but the ratio is proportional to the size of the 
gap, so S" ^ C/T + S as r (similar results have 
been found when applying scaling theory to the Ander- 
son transitior£3) . Here this is "allowed" thermodynami- 
cally, because the "ground-state" of the insulating phase 
has nonzero entropy, since we forced the system into the 
paramagnetic insulating phase. In a real system, how- 
ever, there must be a transition to a ground state where 
the entropy is quenched. In such a case, one would ex- 
pect to see a large peak in the thermopower at low en- 
ergies, with the thermopower ultimately going to zero as 
T ^ 0. We know of no real correlated insulator that has 
a diverging thermopower as T — > 0. 

The computation at very small values of T is diffi- 
cult because the electrical conductivity (or equivalently 
Lii) approaches zero and there are numerical difficul- 
ties associated with properly calculating the conductivity 
in this regime due to the cancellation of two large and 
nearly equal numbers. Since the thermopower requires 
the electrical conductivity, if that cannot be calculated 
accurately then spurious behavior will be seen in the ther- 
mopower. The thermal conductivity looks similar to the 
weaker correlated case, except it appears to go to zero at 
a nonzero temperature which is the expected behavior for 
a correlated insulator with a gap and the potential energy 
pieces become increasingly more important (particularly 
at low temperature). 

In both of these half- filled cases, that of a correlated 



metal and a correlated insulator, the low-temperature 
thermopower is determined by a slightly larger contribu- 
tion from the kinetic energy piece of the heat current than 
from the potential energy piece of the heat current. The 
thermal conductivity, on the other hand has an evolution 
of going from a result nearly completely determined by 
the kinetic energy only piece of the heat current correla- 
tion functions to one where the potential energy pieces of 
the heat current contribute progressively more and more 
to the total thermal conductivity. 

Next we present results for the case where the total 
fiUing Pf. + (w) = 1.5 is a constant but the electrons 
can change from localized to itinerant (i.e., we fix the 
total electron concentration not the individual electrou. 
concentrations). We choose values of the parametersliJ 
where the system has a sharp transition from a state at 
high temperature that has large f-occupancy {{w) w 0.36 
for 0.2 < T < 0.8), to a state at low temperature with no 
f-clectrons (the crossover occurs near T — 0.04). We find 
that the results do not depend too strongly on the param- 
eters in this regime, and choose Ef = —0.7 and C/ = 4 
as a canonical system that is similar to YbInCu4. The 
main difference from the symmetric case studied above is 
that the localized electron filling goes to zero as T — > 0. 
Hence both the kinetic and potential contributions to 
the thermopowcr become small in this limit, and there 
is no large cancellation between two nearly equal num- 
bers to determine the thermopower. We see that in the 
thermal conductivity, the contributions from the mixed 
kinetic and potential energy pieces are the most impor- 
tant, which is an indication of the strengthening of the 
correlations in the system. This feature is hard to see 
from the shape of the thermal conductivity itself. 




O 



§ -1 f-' s(T) - 

^ i' S^T) 

K -2 - s''(T) 

I I I I 

0.5 1 1.5 2 2.5 

Temperature T [t*] 
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FIG. 6. Thermopower for the case U = 4, wi + pd ~ 1-5, 
and Ef ~ —0.7. The main panel shows the two different 
contributions from the kinetic and potential energy pieces of 
the heat current and the inset shows the net thermopower. 
Note the absolute scale for the thermopower is much larger 
here. 




Temperature T [t*] 

Freericks, and Zlatic, Phys. Rev. B, Figure 7 

FIG. 7. Thermal conductivity for the case (7 = 4, 
wi + Pd = 1.5, and Ef ~ —0.7. The plot shows the different 
contributions from the kinetic and potential energy pieces of 
the heat current. 

The thermopower has an interesting exchange of im- 
portance of the different pieces of the heat current as 
temperature is varied. For high temperatures, the ther- 
mopower is determined by both parts of the heat current 
and a compensation effect is important. In the moderate 
temperature regime, the thermopower is dominated by 
kinetic-energy pieces, which then give way to the poten- 
tial energy domination at low temperature, that eventu- 
ally shrinks as T ^ and the thermopower vanishes. 

Note, would need to reverse the sign of the ther- 
mopower to describe YbInCu4, since its charge carri- 
ers are holes rather than electrons. We should also re- 
mark that the thermopower of the Falicov-Kimball model 
does not have a low-energy peak associated with the 
"valence-change" transition — such a sharp peak occurs in 
YbInCu4-like systems because of hybridization effects not 
included in this model. Consequently, the thermopower 
also does not have a low-tcntpfirature sign change as 
seen in the experimental datat3 (a sign change occurs 
at T « 2), but we see that if one could reduce the poten- 
tial energy piece of the thermal current, then the kinetic 
energy contributions to the thermopower could cause a 
sign change to occur. This cannot happen in a pure 
Falicov-Kimball model though, because of the Jonson- 
Mahan theorem. The sign change for the pure Falicov- 
Kimball model generically occurs at much larger values 
of temperature, on the order of the bandwidth. As re- 



gards the thermal conductivity, we find k is dominated by 
the potential-energy piece as T ^ and the mixed piece 
yields a negative contribution over a wide temperature 
range. 

IV. CONCLUSIONS 

We have examined thermal transport in the spin-one- 
half Falicov-Kimball model. We chose this model because 
the transport properties can be solved exactly, and they 
provide an alternate proof of the Jonson-Mahan theo- 
rem for the thermopower. We provide the proof in two 
different ways. The first is a brute-force application of 
the dynamical mean field theory to calculate all relevant 
correlation functions and combine all terms to yield the 
final expressions for the thermal transport coefficients. 
The second is based largely on the techniques of Jon- 
son and Mahan, but one can determine the important 
"generalized polarization" functions exactly in the large 
dimensional limit. Here we extend the Jonson-Mahan ar- 
guments to show analogous results hold for the thermal 
conductivity as well. 

Our formulation also allows us to decompose the con- 
tributions to the thermopower and the thermal conduc- 
tivity into the respective contributions from the kinetic 
energy piece and the potential energy piece of the thermal 
current. We find that generically, these pieces are large 
and opposite in sign for the thermopower so that thermal 
transport carried by the kinetic heat current is almost 
completely compensated by the potential heat current 
producing a small net thermopower. For the thermal con- 
ductivity, we see an evolution of the transport being dom- 
inated first by kinetic energy terms and then potential 
energy terms as the strength of the correlations increase. 
We note, that because the kinetic-energy contribution to 
the thermopower can be straightforwardly determined for 
a number of models, any Hamiltonian that satisfies the 
Jonson-Mahan theorem can be separated into its kinetic 
and potential pieces for the thermopower by simply sub- 
tracting the kinetic energy piece from the Jonson-Mahan 
result. 
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